How loopy is the quantum bounce? — A heuristic analysis 
of higher order holonomy corrections in LQC 
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A well-motivated extension of higher order holonomy corrections in loop quantum cosmology 
(LQC) for the fc = Friedmann-Robertson- Walker model is investigated at the level of heuristic 
effective dynamics. It reveals that the quantum bounce is generic, regardless of the order of correc- 
tions, and the matter density remains finite, bounded from above by an upper bound in the regime 
of the Planckian density, even if all orders of corrections are included. This observation provides 
further evidence that the quantum bounce is essentially a consequence of the loopy nature (i.e. in- 
trinsic discreteness) of LQC and LQC is fundamentally different from the Wheeler-DeWitt theory; 
it also encourages one to construct the quantum theory of LQC with the higher order holonomy 
corrections, which might be understood as related to the higher j representations in the Hamiltonian 
operator of loop quantum gravity. 

PACS numbers: 98.80.Qc, 04.60.Pp, 03.65.Sq 



I. INTRODUCTION 

Over the last years, the status of loop quantum cosmol- 
ogy (LQC) has progressed significantly and has become 
an active area of research. Specifically, with the inclusion 
of a free massless scalar field, the comprehensive formu- 
lation of LQC in the k = Friedmann-Robertson- Walker 
(FRW) (i.e. spatially flat and isotropic) model has been 
constructed in detail, giving a solid foundation for the 
quantum theory and revealing that the big bang singu- 
larity is resolved and replaced by the quantum bounce, 
which bridges the present universe with a preexisting one 
0, 0, S] ■ Resolution of the classical singularity and oc- 
currence of the quantum bounce have been shown to be 
robust and similar results are also affirmed for ex- 
tended models [IlEHIi. 

However, despite an attractive and long-sought fea- 
ture, it is questionable whether the quantum bounce re- 
sults intimately from the quantum nature of Riemannian 
geometry of loop quantum gravity (LQC), as the same 
result can be easily obtained even at the level of heuris- 
tic effective dynamics without invoking the sophisticated 
features of LQC [1, 0, HH- 

To quickly see this, we start with the classical Hamil- 
tonian constraint for the A; = FRW model: ^ 



H = H,r 



(1) 



where c and p are the Ashtekar variables satisfying the 
canonical relation 



{c,p} = 



87rG'7 



(2) 
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^ More precisely, p should be \p\ in JTJ. The sign of p corresponds 
to spatial orientation, which is irrelevant for our purpose. 



with 7 being the Barbero-Immirzi parameter and the 
conjugate momentum of the free massless scalar field 

4>{x, t) — (j){t) with 



1. 



(3) 



Next, at the heuristic level, we take the prescription of 
"holonomization" to replace c with 



sm{fic) 



(4) 



by introducing the discreteness variable p,. The heuristic 
effective dynamics is then solved as if the dynamics was 
classical but governed by the new "holonomizcd" Hamil- 
tonian, which reads as 



Hp 



sm fic 



2p3/2 • 



(5) 



Particularly, in the improved dynamics suggested by 
/i is given by 



(6) 



where A is the area gap in the full theory of LQC and 
A = 2'\/37r7£p[ for the standard choice (but other choices 
are also possible) with £pi :— \J Gh being the Planck 
length. With ^ imposed, the modified Hamiltonian con- 
straint Hf^ — immediately sets an upper bound for the 
matter density: 



P<P ■= 



■ sin < 3p 



2p^ SirGjA ' 
where the Planckian density is defined as 
PPI :- {STTG^^AyK 



PI, 



(7) 



(8) 



Apparently, without going into the detailed construc- 
tion of LQC at all, it is anticipated that the matter den- 
sity is bounded above and thus the quantum bounce is 



2 



expected.^ One might then argue that the boundedness 
of has httle to do with the fundamental structure of 
LQC but merely results from the formal modification of 
In this sense, the quantum bounce seems to be an 
ad hoc phenomenon and not really "loopy" enough. 

In response to this criticism, a simplified but exactly 
soluble model of LQC has been studied and used to show 
that the quantum bounce is generic, not just restricted to 
the states which are semiclassical at late times Fur- 
thermore, the study of [J] brings out the precise sense 
in which the Wheeler-DcWitt (WDW) theory approxi- 
mates LQC and the sense in which this approximation 
fails, thereby showing that LQC is intrinsically discrete 
and the underlying discreteness is essential for the quan- 
tum bounce. 

To add further evidence for the loopy nature of the 
quantum bounce, we explore a new avenue by investigat- 
ing the well-motivated extension of higher order holon- 
omy corrections, with which the prescription of holon- 
omization is more involved than ([4]), yet the quantum 
bounce is still ensured, at least at the level of heuristic 
effective dynamics. 



II. HIGHER ORDER HOLONOMY 
CORRECTIONS 

One of the very features of LQC is that the connec- 
tion variable c does not exist and should be replaced by 
holonomies. Following the standard techniques in gauge 
theories, components of the curvature F — TkPt'jdx'^ Adx^ 
can be expressed in terms of holonomies (i.e. Wilson 
loops) as 



rTr 



i^l - >) 



(9) 



where 2iTk = <Jk are the Pauli matrices and hi^^ is the 
holonomy around the square whose edges are par- 
allel to the i and j directions and of coordinate length 
p,L {L is the coordinate length of the edges of the fidu- 
cial cell V prescribed to make sense of the Hamiltonian) . 
This amounts to approximating in ([1]) by sin^ pc/ 
and thus the prescription of (For more details, see 
Appendix B of [ll|-) 

In the standard treatment of LQC, instead of shrinking 
to a point, the expression ^ is regarded as funda- 
mental and fjL is set to a nonzero value as an imprint of 
the discrete area spectrum in the full theory of LQC. 
It is natural to ask whether the approximation ^ can 
be improved; that is, is it possible to approximate the 



^ If a different scheme otlier than is adopted, can also 
be shown to bounded above, as long as p, is prescribed to be 
p, oc {A/pY with r > — 1, but the upper bound depends on the 
constant of motion p^ except the case of r = 1/2. 



curvature F in terms of holonomies along edges of finite 
lengths, yet with arbitrary accuracy? 

The approximation ([9]) is based on the Stokes' theorem 
and, in general, it becomes exact only in the limit when 
Dij shrinks to a point. In the context of cosmologies, 
however, thanks to homogeneity, for any given finite fl, 
it is still possible to approximate c in terms of sin fic to 
arbitrary accuracy. As a heuristic approach, disregarding 
the standard Stokes' theorem but instead considering the 
Taylor series 



sin ^ X — ^ 



{2ky. 



k=0 



22'=(fc!)2(2fc-f 1 



■ X 



2k+l 



(10) 



for — 1 < X < 1 and setting x — sin /ic, we have 

1 ^ (2fc)! . ,2k+i ,,,, 

This inspires us to define the nth order holonomized con- 
nection variable as 



mi 



1 " 

22'=(fc!)2(2fc + 1) 



(sin fic) 



2k+l 



(12) 



which can be made arbitrarily close to c (as n — )■ oo) but 
remains a function of the holonomy s'm flc and the dis- 
creteness variable p.. Therefore, to implement the under- 
lying structure of LQC by replacing c with holonomies, 
cjj"' can be used as an improved version of ^ , which now 

reads as c|j"~°^ . 

Instead of (O , the Hamiltonian with holonomy correc- 
tions up to the nth order can be designated as 



87rG7 



2p3/2 ■ 



(13) 



If we take ([131) the departing point for LQC, the quan- 
tum theory is much more difhcult to construct, but it 
might be possible to treat the higher order corrections as 
perturbations based on the well-established n = formal- 
ism. The higher order corrections correspond to higher 
powers of sin /Ic, which might be understood as the im- 
print of generic j representations for holonomies in the 
Hamiltonian operator in the full theory of LQG [T^. [T3|. 
(See also [l3| for the issues of j ambiguity in LQC.) Thus, 

(n) 

even though cj^ is obtained heuristically, the modified 
Hamiltonian of (jl3p may reflect the underlying physics 
of LQG in a more elaborate fashion. 

One might suspect that if we include corrections of all 
orders, the quantum theory of H^~°°^ will lead to the 

same result of the WDW theory, as c^""-* — c formally. 
This should not be the case, because the elementary vari- 
ables are still sin pc and p, instead of c and p, even in the 
limit n — > oo and, therefore, the striking difference be- 
tween LQC and the WDW theory as emphasized in [3| 
should persist. 
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While the quantum theory of LQC could be extremely 
difficult, we will investigate the ramifications of the 
higher order holonomy corrections at the level of heuristic 
effective dynamics. Beforehand, using we compute 



{c,ci")} 
feci")} 



87rG7 dfi 
3/2 dp 
8nGj 



cos(/ic) 6„(/ic) c 



cos(/ic) 6„(/ic), 



where 



6„(^c) 



E 

fc=0 



|cos /ic| 



(2fc)! , . _ , 



2fc 



(14) 
(15) 

(16) 
(17) 



Remark. One should not confuse the extension of 
higher order holonomy corrections discussed here with 
that studied in [1^ and [l^. The former extends 2^""°^ 

(n) .... . _ . 

to Cf^ with the inclusion of higher powers of sin fic, which 
are most likely to be interpreted as higher j representa- 
tions, while the latter takes into account the error cor- 
rections in © with respect to the powers of p,. That 
is, the former is about the O (sin^""*" ^ /ic) corrections for 
c while the latter is about the ©(/i^"^^) corrections on 
sin^ fic. These two extensions are motivated differently 
and give distinct dynamical behaviors. As will be seen, at 
least for heuristic effective dynamics, the study of this pa- 
per afhrms that, for any order n, the nonsingular bounc- 
ing scenario is generic, as opposed to the result of [T5| 
and [III, which indicates that the 0{p,'^) corrections lead 
to the nonperturbative effects (hyper-inflation / deflation) 
and thus yield a qualitatively different cosmological sce- 
nario in which the bounce and singularity coexist. More 
comments on the distinction are given in Sec. IIVI 



III. HEURISTIC EFFECTIVE DYNAMICS 

At the level of heuristic effective dynamics, the evolu- 
tion is solved as if the dynamics was classical but gov- 
erned by the new Hamiltonian (|13p with holonomy cor- 
rections up to the nth order. For the case of n = 0, it 
has been shown that this heuristic treatment gives a very 
good approximation for the quantum evolution of LQC, 
and the bouncing scenario of the effective solution gives 



the absolute upper bound for the matter density [1 



For n > 0, the reliability remains to be justified, but the 
heuristic analysis can still provide good ideas of what 
the quantum evolution may look like in the presence of 
higher order holonomy corrections. 

By choosing the lapse function N = p^^^ associated 
with the new time variable t' via dt' = N^^dt [t is the 
proper time), the modified Hamiltonian (|13p is rescaled 
and simplified as 



(»)' 



(4' 



(»)• 



2 2 
P 



(18) 



The effective equations of motion are then given by 
Hamilton's equations: 



dt' 

d(l) 

dt 
dc 

dt' 



dp 
dt' 



{p^,Hjr^'}^o ^ p, 

r iir(")'i ^ (")^2 



is constant, (19) 
(20) 
(21) 



2 1 9/2 

J p, dp 



(n) 



cos(/ic) 6„(/ic) c — c)j 



{p, H';" } - - cos(/2c) 6„(/2c) 4")p2^ (22) 
7 



(") 2 
4 P , 



and the constraint that the Hamiltonian must vanish: 



pI 



2 87rG7 
In particular, these lead to 



—,ict'W- (23) 



1 dp 

p dcj) 



IQttG 



cos(^c) 6„(/ic), 



(24) 



as (f) is treated as the internal time. 

In the classical regime, we have /ic ^ 1 and thus 
cos(/ic) 1, sin(/2c) 0, (3n(/ic) — > 1 and cj^^ c; 
therefore, the above equations all reduce to their classical 
counterparts. In the backward evolution, the quantum 
corrections are more and more significant as /2c becomes 
appreciable. Eventually, p gets bounced at the epoch 
when cos(/2c) in ([M)) flips signs. The exact point of the 
quantum bounce is given by cos(/2c) = (i.e. /2c = 7r/2). 
By ([T2|) . this happens when 



(2fc)! 



fc=0 



22fc(fc!)2(2/s + l) 



— ■ ■Sn 



(25) 



(note that ^ 7r/2 as n — > oo), or equivalently, by (|23p. 
when 

P.p = Pcl--=^SlPPi (26) 
if the improved scheme ([6]) is adopted. Note that 

37r2 



3pp, = p^^l < f^, <■■■< p)-> = :^pp, (27) 



This shows that the occurrence of the quantum bounce 
and the boundedness of the matter density are generic for 
any given n. Even at the limit n — > oo, the matter density 
remains bounded, as opposed to the classical theory. 
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100 I 1 1 1 ^ 1 . 1 . 1 1 I 1 1 1 1 1 1 1 1 1 

-2.0 -1.8 -1.6 -1.4 -1.2 -2.0 -1.8 -1.6 -1.4 -1.2 

FIG. 1: Solutions of the heuristic effective dynamics with = 2. x 10''ft\/87rG, the initial value p{(t>o) = 1. x 10^£pi and 
7 = ln2/(\/37r). (a) Solid curves represent p{4>) for n = 0,1,2,3,4,10,20,30,40,50 (the larger n is, the smaller p is at the 
bouncing point). The dotted line is the classical solution, (b) Solid curves are the corresponding p^, the peaks of which agree 
with p^";' given by (|26p (the larger n is, the bigger p'";). is). The solution of n = oo is also indicated as the dotted curve, which 
is sharply kinked at the top tip at the value of p'^'^^ given by (|27p . 



IV. DISCUSSION 

For a given initial condition which satisfies the Hamil- 
tonian constraint l|23|) . the differential equations (|19p- 
([^ can be solved numerically (by the Runge-Kutta 
method) to yield the detailed evolution of the heuristic 
dynamics for a given n. The numerical solutions with 
^ imposed are depicted in Fig. [l]for different n.^ It is 
affirmed that the nonsingular bouncing scenario is robust 
regardless of n: Two classical solutions (expanding and 
contracting) are bridged by the quantum bounce, which 

takes place when approaches the critical value 

As n increases, the upper bound p^"j{ increases as well, 

but the limiting value p^^j'* remains finite in the order of 
the Planckian density as indicated in P7)) . Furthermore, 
in Fig. [T]^a) it is noted that the larger n is, the longer 
the effective solution follows the classical trajectory be- 
fore being deviated by the bounce. Put differently, the 
quantum bounce takes place more abruptly for a larger 
n. In the extreme case of n oo, the effective solution 
exactly matches the classical solutions on both sides of 
the bounce, which is so abrupt that it only imprints a 
kink on the solution of p{(t>). This is expected, since for- 
mally c|j"'* approximates c closer and closer as n increases, 
but nevertheless the factor on the right-hand side of (IM)) 



The double precision used in the computer program loses neces- 
sary accuracy and results in noticeable round-off error when n is 
huge (n > 100); thus, the numerical computations are done only 
up to n = 50. 



yields the limit: cos(/2c) S„(/2c) cos(/lc)|cos(/2c)|~ — 
sgn(cos(/ic)), which gives rise to the quantum bounce as 
a kink. The remarkable point is that, even in the limit 
n — » oo, the heuristic effective dynamics does not simply 
reduce to the classical one. 

If the quantum theory of LQC with holonomy correc- 
tions up to the nth order can be constructed, following 
the lesson of 3, we expect that the expectation value 

of the matter density remains finite as p^"j{ obtained in 
([26|) sets the absolute upper bound. Even if all orders of 
holonomy corrections are included, it is very likely that 
the quantum bounce persists and the expectation value 
of the matter density remains lower than p^^^.'' given by 
(|27p . This suggests that, even with n — s- oo, the quantum 
theory of LQC, if it can be constructed, is fundamentally 
different from the WDW theory. Furthermore, in the 
quantum theory of LQC of n — > oo, the abruptness of 
the bounce could be well smoothed by the quantum fluc- 
tuations. 

As higher orders of holonomy corrections are included, 

(n) 

c\ gets closer to formally agreeing with c, yet the 
heuristic analysis shows that the matter density remains 
bounded and the bouncing scenario holds for any arbi- 
trary n. This makes it more convincing to assert that the 
occurrence of the quantum bounce is essentially a con- 
sequence of the loopy nature (i.e. intrinsic discreteness) 
of LQC. The fact that it is well motivated and makes 
sense to incorporate higher order holonomy corrections 
at the level of heuristic effective dynamics encourages one 
to construct the quantum theory of LQC based on the 
new Hamiltonian (|13p . the investigation of which might 
in turn shed light on the issues of j ambiguity for the 
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Hamiltonian operator in the full theory of LQG. 

The observation that the bouncing scenario survives 
the higher order holonomy corrections is also expected 
for other extended models. With the inclusion of generic 
matters, for example, the resulting effective solutions 
should be qualitatively the same as those of the familiar 
n = formalism, since the matter part of the Hamilto- 
nian constraint involves only p but no c and hence re- 
ceives no holonomy corrections. Consequently, the sce- 
nario in Appendix A of [Slwith a nonzero cosmologi- 
cal constant and those in p, Q for the k — ±1 FRW 
models should live on, because at the level of effective 
dynamics the cosmological constant and the spatial cur- 
vature terms can be treated as matter with the equation 
of state parameters w = —1 and —1/3, respectively. In 
the context of anisotropic models, the same is also antic- 
ipated, but the exact value of the absolute upper bound 
for the matter density is difficult to pinpoint, as the ef- 
fective equations of motion are more complicated and 
the exact condition for the bounce depends on the de- 
gree of anisotropy (see [ll| for the Bianchi I model and 
for the Kantowski-Sachs spacetime for the case of 
n = 0). Similarly, the loop quantum geometry of the 
Schwarzschild black hole interior predicted in [20| is not 
expected to be spoiled either, but again the details de- 



mand closer examination. 

Additionally, as remarked earlier in Sec. [Ill it should 
be noted that the extension of higher order holonomy 
corrections studied in this paper is distinct from that in 
[l^[3l. Opposed to our conclusion, [15] and [3| sug- 
gest that the nonsingular bouncing scenario that appears 
in the lowest order could only be an artifact of simplifi- 
cation, since it is qualitatively modified by the 0{fl'^) 
corrections. As it is debatable which extension approach 
makes more sense, whether the higher order holonomy 
corrections modify the bouncing scenario qualitatively or 
only quantitatively remains an open question. To give a 
decisive answer, construction of the quantum theory of 
LQC with the inclusion of higher order holonomy correc- 
tions and further investigations from the perspective of 
LQG are necessary. 
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